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Abstract 

Continuous phase transitions are studied in a two dimensional nonequi- 
librium model with an infinite number of absorbing configurations. 
Spreading from a localized source is characterized by nonuniversal crit- 
ical exponents, which vary continuously with the density (ft in the sur- 
rounding region. The exponent 5 changes by more than an order of 
magnitude, and rj changes sign. The location of the critical point also 
depends on (ft, which has important implications for scaling. As expected 
on the basis of universality, the static critical behavior belongs to the 
directed percolation class. 
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I. INTRODUCTION 



Despite recent progress, understanding the phase diagrams and critical points 
of nonequilibrium models remains a challenging problem. In the relatively well- 
understood case of continuous transitions into an absorbing state (i.e., one that traps 
the dynamics), a high degree of universality has been found, supporting the pre- 
diction that such transitions belong generically to the class of directed percolation 
(DP) [|T]-[J. Examples include the basic contact process and its variants sur- 
face reaction models P JiTl|JTT[| , branching annihilating random walks with odd parity 
121-141, and assorted multiparticle processes [1^-|19|. The static critical behavior of 



models with multiple absorbing configurations PQ -|27| also falls in the DP class. On 
the other hand, some aspects of time- dependent critical behavior in the latter class 
of models are surprising: the critical exponents which describe spreading from a seed 
vary continuously with the density in the environment, and the usual hyperscaling 
relation must be modified. Until now this phenomenon has been investigated in one 
dimension only; here I present the first study of nonuniversal critical spreading in a 
two dimensional model. 

Models with an absorbing state are characterized by an order parameter p — the 
density of self-reproducing entities or active sites. The stationary value p vanishes as 
the reproduction rate A approaches a critical value A c : p oc A 13 for small A = A — A c . 
For A < A c the only steady state is the absorbing state, in which all change ceases. 



Much insight is gained from studies of spreading from a single active site p8| . The 
chief quantities of interest are the survival probability, P{t), (the process is said 
to survive so long as it has not become trapped in the absorbing state), the mean 
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number of active sites n(t), and the mean-square distance, R 2 (t), of active sites from 
the original seed. At the critical point they follow asymptotic power laws, 

p(t) oc r s , (l) 



n{t) oc f, (2) 

and 

R 2 (t) oc t z . (3) 

These power laws, which form the most salient feature of the critical spreading process, 
have been confirmed for a broad class of models. The exponents 6, rj, and z, estimated 



from simulations and series expansions [2E] are in good agreement with DP values. 



Several models have been devised with absorbing states that embrace a large num- 



ber of configurations: two-dimensional catalysis models p2|-p5|, and one- dimensional 
particle systems such as the pair contact process (PCP) p0|j21| , |27|] . In the PCP each 
site of a lattice (Z) is either vacant or occupied by a particle; the active "sites" are 
nearest-neighbor particle pairs. They are annihilated with probability p, and give 
birth, with probability 1 — p, to a new particle at a vacant neighboring site. Since 
particles may only be created or destroyed if there are pairs, and since particles cannot 
move to other sites, any configuration of isolated particles is absorbing; the system 
always becomes trapped in such a configuration if p > p c . While static critical be- 



havior in models with multiple absorbing configurations still conforms to DP [[30 
spreading presents a new feature: one can choose the absorbing configuration into 
which the seed is placed. Critical spreading has been studied in three rather dif- 
ferent one-dimensional models with multiple absorbing configurations: the PCP, the 



threshold transfer process, and the dimer reaction |21|j27[] . For "natural" absorbing 
configurations (those generated by the critical process itself), the spreading exponents 
assume DP values. If, however, the initial density 0; differs from the natural value, 
the exponents 5 and rj change; they vary linearly with fa. (<j) represents the density 
of isolated particles in the PCP, not the order parameter.) A new exponent, f3' must 
be introduced to describe the ultimate survival probability: 

Poo = lim P(t) ex (p c - pf (4) 

in the active regime. This exponent also depends on fa. (j3' = j3 for natural initial 
configurations, just as for models with a unique absorbing configuration.) It appears 
in the generalized hyperscaling relation [27|, 



:i + f)^ = y, (5) 



which reduces to the relation originally derived by Grassberger and de la Torre |28| 
for f3' = (3. Certain properties are independent of fa: the location of the critical point, 
the exponent z, and 5 + rj, which governs population growth in surviving trials. This 
suggests that the initial density influences the survival probability, reflected in the 
exponents 5 and f3\ but does not alter the asymptotic properties of surviving trials. 

The background sketched above motivates the investigation of critical spreading 
in dimensions > 2, which offers a greater range of possible growth patterns. In the 
following section I define a two-dimensional model with multiple absorbing configura- 
tions, and describe the simulation procedure used to study it. Results are presented 
in Sec. Ill, and a scaling analysis in Sec. IV, followed by a brief summary in Sec. V. 
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II. THE MODEL 



The second neighbor reaction (SNR) is a Markov process or interacting particle 
system in which each site of the square lattice Z 2 is either vacant or occupied by a 
particle. Particles may not occupy adjacent sites. If all four nearest neighbors of a 
vacant site (i,j) are also vacant (i,j) is said to be open. Open sites are the only active 
sites in the model. In each step of the process, a site is chosen at random; if it is 
open, a particle is placed there provisionally, otherwise nothing happens. Suppose a 
particle has just arrived at (i,j). If all four of the second neighbors (i±l, j±l) are 
vacant the new particle is adsorbed, and remains at (i,j). But if any of the second 
neighbors are occupied, the new particle is only adsorbed with probability p; with 
probability 1 — p it reacts with one of its neighbors (selected at random, if necessary), 
and both particles leave the lattice. Each such reaction creates one or more new open 
sites. Adsorbed particles are immobile, and can only leave the surface in a reaction 
with a newly-arrived particle. Each adsorption attempt marks a fixed time interval; 
we take the time unit as N = L 2 such events. 

Since the fraction of open sites is quite small in the vicinity of the critical point, 
efficiency is greatly improved by choosing trial adsorption sites from a list of open 
sites. This entails a variable time increment, At = 1/N a , for each event, where N a is 
the number of open sites. Taking statistics at fixed time intervals (rather than after 
a fixed number of events), ensures proper stationary averages. 

For p — 1 the model corresponds to random sequential adsorption (RSA) with 
nearest-neighbor exclusion; starting from an empty lattice, particles adsorb until no 
open sites remain. In this case there is an exponential approach to the jamming 
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density 0oo = 0.36413 f31| , |32|| . The order parameter p for the SNR is the fraction 
of open sites; configurations devoid of open sites are absorbing. Of course, many 
absorbing configurations are possible, with particle densities (f) ranging from 1/5 to 
1/2 (see Fig. 5). I shall refer to the set of absorbing configurations characterized by 
common statistical properties (particle density and correlations) as an absorbing state. 
It is also useful to define the ensemble A Pj l of absorbing configurations generated by 
the SNR with adsorption probability p, starting from an empty lattice of size L. For 
any absorbing configuration £ on L 2 sites, A p ,l{() ls the probability of realizing ( as 
the final configuration. 

Based upon experience with models of this kind — in particular, with the dimer 
reaction, the one-dimensional cousin of the SNR |2l].[27]] - one expects that there 
is a critical value p c above which the system always gets trapped, whilst for p < p c 
there is an active state in which the reaction may proceed indefinitely. Equivalently, 
one expects that for p > p c the absorbing state is the only stationary state, and that 
for p < p c there is in addition a nontrivial invariant measure on configuration space. 
(While the only latter emerges in the infinite-size limit, for p < p c the quasi- stationary 
state is sufficiently long-lived to allow precise characterization in simulations.) 



III. SIMULATION RESULTS 



An accurate estimate of p c is crucial for reliable determination of critical behavior. 
Several simulation methods are available for locating the critical point: half-life stud- 
ies, time-dependent simulations, and analysis of stationary properties. Use of any one 
method exclusively can be misleading. Here I begin with the half-life method, which 
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is based on finite-size scaling analysis of the mean survival time. One generates a 
large number of independent realizations, all starting from an empty lattice (all sites 
open), with periodic boundaries, and determines the time r required for half the sam- 
ple to reach the absorbing state, r is a decreasing function of p and an increasing 



function of lattice size L. One expects that in the vicinity of the critical point |TS] 



r(p, L) ~ L v ^ v± h(AL 1 / ]/J -) , (6) 

where A = p c — p and h is a scaling function. The exponents v\\ and u± govern 
the divergence of the relaxation time t r and the correlation length £ as A — > 0: 
t r oc A~^n and £ oc A~ UJ -. The critical point is distinguished by the simple power- 
law r oc L v \\l Vl - as L — > oo. Following a preliminary survey which indicated that 
p c ~ 0.3935, I estimated r for L = 16, 32, 48 64, 96, 128, and 176, for p = 0.3930 
- 0.3940. It is convenient to plot lnr — (v\\/v±) hxL versus InL, using the expected 
value of the exponent ratio for DP in 2+1 dimensions (see Fig. 1). Only for p in 
the range 0.3926-0.3927 are the data consistent with a power law. The corresponding 
exponent is uu/u± = 1.76(2), in good agreement with the DP value, 1.764(26). 

More precise estimates of p c may be obtained from simulations of spreading |[28| . 
As noted above, the critical point is marked by asymptotic power laws, Eqs. (|IJ) - (Q). 
This part of the study focuses on the natural absorbing state — configurations drawn 
from A Pci l, defined in Sec. II. Initial configurations are generated by running the 
simulation (starting with all sites vacant), at a given p, on a lattice of size V = 128, 



until it reaches an absorbing configuration. Following the procedure of Ref. [^T 
copies of this configuration are assembled into an absorbing configuration of size 
L = 1280. Then a single (randomly chosen) particle is removed, creating an open site 



(or possibly a small cluster of open sites), which serves as the seed for the spreading 
process. 

In spreading simulations the process halts before reaching the boundary; there are 
no finite-size effects as such, though one does of course expect finite-time corrections 
to scaling. Trials ran to a maximum time t max = 10 4 ; 5 x 10 4 independent trials were 
generated for each p value studied. The survival probability, mean number of open 
sites, and mean-square spread for p = 0.392615 are plotted in Fig. 2, showing clear 
evidence for power laws. In order to fix p c and the exponents with precision, I plot, 
in Fig. 3, the local slopes of the log-log graphs versus 1/t, for several values of p. 
(The local slope St is found from a linear least-squares fit to the data for lnP(t) over 
the range [In t — 1.5, \nt + 1.5], and similarly for the other exponents.) I conclude 
that 0.39261 < p c < 0.39263, given the marked curvatures of the St and r] t graphs for 
p-values outside this range. The exponents obtained by extrapolating the local slopes 
to infinite-t are in good agreement with accepted values for DP in 2+1 dimensions 
(see Table I). 

I turn to the exponent j3' governing the ultimate survival probability starting from 
a localized active region. Natural absorbing configurations, constructed as above 
(with p = p c = 0.39262), were used in determinations of for p < p c . These studies 
yield (3' = 0.58(2), whilst in DP (3 = (3' = 0.586(14) |33J. The stationary open site 
density p was determined on lattices of size L = 16, 32, 64, and 128. To minimize 
finite size effects, p(p) was obtained for two or more system sizes, until doubling L 
produced no significant change. A linear least-squares fit to a plot of In p versus 
ln(p c — p), (with p c = 0.39262), yields f3 = 0.56(2). In summary, the static critical 
behavior and critical spreading from "natural" absorbing configurations are both fully 



8 



consistent with DP scaling in 2+1 dimensions. 

It is of interest to characterize the natural absorbing configurations. Large sam- 
ples of the latter were generated by running the simulation at p c , starting from an 
open lattice (L < 176), until it reached an absorbing configuration. The results for 
the particle density cj) na t show only a weak dependence on L, permitting reliable ex- 
trapolation to infinite-L, yielding cj) nat = 0.33429(3). Further insight into absorbing 
configurations is provided by the two-point total correlation function, 

h(r)5 < g Mff + r » -l, (7) 

where cr(x) is an indicator variable which is 1 iff site x is occupied. Fig. 4 shows that 
h(r) is negligible for r > 6. Thus correlations in natural absorbing configurations in 
the SNR are of only slightly greater range than in the one- dimensional dimer reaction 

Having established the properties of the active stationary state, the natural ab- 
sorbing state, and of critical spreading from the latter, I turn to spreading in "non- 
natural" absorbing configurations. Detailed studies were performed for the maximal 
density {(f) = 0.5), the minimal one (0 = 0.2), and for two kinds of configurations 
with <p = 0.25, one with square-lattice symmetry, the other having a preferred axis. 
In addition to these periodic states, absorbing configurations were generated via RSA. 
(Examples are shown in Fig 5.) In each case, the critical point, p c {4>), was identified 
by searching for asymptotic power-laws as in Eqs ([[]) — (|3]). 

Consider first the lowest density absorbing configuration (<f) = 0.2), in which 
particles are arrayed in a y/E x \/5 pattern. Simulations indicate that the spreading 
process is supercritical for p = p c ,buik = 0.39262. The survival probability P(t), for 
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example, tends to a nonzero value as t — > oo, rather than exhibiting power-law decay. 
Power-law spreading occurs for p = 0.4102(2), about 4% greater than p c ,buik (see Fig 
6). As indicated in Table I, the spreading exponents rj and z are much larger than the 
corresponding DP values, and 5 and j3' are much smaller. The changes in p c and in 
the exponents accord with the notion that open sites spread more readily in a region 
of low particle density. 
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TABLES 

TABLE I. Critical point and exponents in the second neighbor reaction. Figures in 
parentheses denote uncertainties. 



1 


Pc 


c 




V 


z 




0.2 


0.4102(2) 


0.078(8) 


0.64(4) 


1.84(2) 


0.125(5) 


0.25 (isotropic) 


0.3965(1) 


0.054(2) 


0.90(1) 


1.76(1) 


0.18(5) 


0.25 (aniso.) 


0.39238(2) 


0.677(5) 


-0.035(5) 


0.97(1) 




0.33443 (natural) 


0.39262(1) 


0.468(6) 


0.216(8) 


1.117(13) 


0.58(1) 


0.36413 (RSA) 


0.39240(5) 


0.64(1) 


0.035(20) 


1.038(3) 


0.70(2) 


0.5 


0.39220(5) 


1.30(5) 


-0.12(6) 


0.91(2) 


2.0(2) 


DP 




0.460(6) 


0.214(8) 


1.134(4) 


0.586(14) 



11 



The absorbing configuration with maximal particle density (0 = 1/2), by contrast, 
presents hostile terrain for spreading. Running at the bulk critical point, the process 
appears to be subcritical: P(t) and n(t) decay exponentially. To observe power-law 
spreading one must reduce p to 0.3922, about 0.14% below the bulk critical value. 
The shift in p c is small, but the effect on the exponents is dramatic: spreading 
is sub-diffusive (z < 1), and rj is negative! (The latter is consistent with known 
constraints: rj + 5 > 0, so the population in surviving trials increases with time.) 
Studies with intermediate initial densities, summarized in Table I and Fig. 7, confirm 
the tendencies noted for the extreme cases: higher densities correspond to smaller 
values of rj and z, larger 5 and f3', and a depressed p c , and conversely, consistent with 
the expectation that the larger 0, the more spreading is impeded. The only exception 
is the anisotropic system, which despite having < (p na t, resembles isotropic systems 
with > (pnat- Why anisotropy renders spreading more difficult is unclear. The 
contrast between the two studies at = 0.25 is nonetheless striking, and demonstrates 
that factors other than the particle density can influence critical spreading. 



IV. SCALING THEORY 



We have seen that as 0, the density of particles in the environment, is varied, 
the location of the critical point changes as well. The shift in p c implies that critical 
processes starting from non-natural absorbing configurations do not evolve to the 
bulk critical state, which exists uniquely at p c ,buik- As a result, the usual scaling 
theory for absorbing states must be modified. In (isotropic) absorbing configurations 
with < (p n at, critical spreading occurs at p = p c > p Ct buik, where an active steady 
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state is not possible. As the population of open sites spreads, it leaves in its wake a 
region devoid of open sites, with particle density 0/, slightly greater than (p nat . After 
some time, this region cannot be reactivated. (For p > p c ^ na t, the survival probability 
decays exponentially even when (ft = cj) nat .) Thus for p c ,nat < V < Pc(<f>), the active 
region is confined to an expanding ring; it is a kind of "chemical wave" which converts 
one type of absorbing configuration to another as it passes. An example of this kind 
of spreading is shown in Fig 8. The mean density of open sites (averaged over 500 
trials), is plotted in Fig. 9 as a function of distance from the origin, confirming that 
the active region forms an expanding ring. Critical spreading with > <fr nat presents 
a rather different picture. Here the critical point lies below p c ,buik, which means that 
surviving trials consist of an expanding region with a small but finite density of open 
sites. 

For non-natural initial conditions, the exponents violate the hyperscaling relation, 
Eq (|5|). This is not surprising, since the scaling argument assumes that surviving 
processes evolve, for large t, into a bulk critical system. [0,0 This assumption, as 



we have seen, is not valid when p c ^ p c ,buik- To describe this novel situation, the 
scaling analysis can be modified as follows. We assume, as usual, that spreading may 



be described in terms of a pair of scaling functions, defined via [[28 

p{x,t) ~t"- <fa / 2 G?(x 2 /t*,At 1 / ,/ ii), (8) 

and 

P(t)^t- s ^{At x ^\), (9) 

where A = p c {4>) — P is the distance from the critical point. (For A > spreading 
may continue indefinitely.) In Eq. p{x,t) is the local order-parameter density, 
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averaged over all realizations; it is concentrated near the origin at time zero. The 
prefactors are constructed to yield power laws — Eqs. ([[]) - © — when A = 0. 
Existence of the limit P^ implies that Q(x) ~ x^' as x — > oo, with f3' = 8v\\. 

Consider spreading with < <p nat and p c (4>) > Pcfiviki m which the active region 
is an expanding ring. Simulations indicate that for A > 0, the active region expands 
at a constant speed v(A), and that the maximum open-site density in surviving 
trials attains a steady value p (A) for large t. In the contact process and similar 
models, p(x, t) —>■ PooP as t — > oo, for any fixed x. But here we expect instead that 
p(x, t) — > P 00 p (A)f(\x.\/vi), where f(u) attains its maximum value (unity) at uq ~ 1, 
and / ^0as«-^0or«-» oo. We expect po to have a power-law dependence on A: 

Po ~ A' 3 ". (10) 



Consider the limit t — > oo, |x| — > oo, with |x|/vt w m . For A small but positive, 
p(x, t) ~ A^'po, so that G(oo,y) ~ yP'+P" for large y. On the other hand, we 
must have G(oo,y) ~ y- u \\(v- dz / 2 ) f or hrn| x i t _ +0O p(x, t) to exist. Comparing these 
asymptotic behaviors, we find a hyperscaling relation for annular spreading 

Inserting the exponent values shown in Table 1 yields j3" = 1.80(35), for cf) = 0.2. 
This is consistent with f3" = 1.68(6), obtained directly from simulations at p- values 
between p C)bu ik and p c (0.2). 

As noted above, surviving trials at p c (4>) for (f> > (p nat develop a finite open-site 
density and so are compact objects. The critical exponents for spreading of compact 



colonies are expected to satisfy the relation pl |: 
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which simply expresses the growth of population in a region of positive density, given 
that the radius grows ~ t z ^ 2 . The exponents for density = 0.5 are, however, 
inconsistent with Eq. (0): 5 + rj — z = 0.27(13) (The deviation from the generalized 
DP hyperscaling relation, Eq. is considerably larger: (l+f3/(3')5+r)—z = 0.6(2).) 
For RSA initial configurations (0 = 0.36413) the corresponding deviations are -0.3(1) 
and 0.2(1). We can understand this discrepancy by noting that the onset of compact 
growth, which requires that the colony diameter D >> £, the bulk correlation length, 
should only occur at very long times, since \p — p c ,buik\/Pc,buik is ver Y small (less than 
10~ 3 ). This is supported by the observation that even at fairly late times (t = 10 4 ), 
the distribution of open sites does not appear compact, but rather consists of several 
disconnected regions (see Fig. 10). Presumably the spreading exponents found for 
= 0.5 and 0.364 are not the asymptotic ones; their determination will require studies 
of considerably longer duration. 



V. SUMMARY 



The first study of nonuniversal spreading in two dimensions reveals a large varia- 
tion of critical exponents with the density of the environment into which the process 
grows. The variation of the spreading exponents is much more pronounced than in 
one-dimensional models with multiple absorbing configurations pl"U2"7|. In the SNR 6 
and f3' vary by more than an order of magnitude, rj changes sign, and z, which is es- 
sentially constant in the one-dimensional models, varies by more than a factor of two. 
(Similarly, 5 + rj appears constant in one dimension, but varies substantially here.) 
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Moreover, the critical point, p c , depends on 0, whereas in one dimension it remains 
constant. The shift in p c has important consequences for scaling since the critical 
spreading process does not evolve into a bulk critical state for <p ^ <t> n at- Compari- 
son of spreading in two environments with the same density but different symmetries 
reveals that the process is strongly influenced by factors other than the density. 
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Figure Captions 

Fig. 1. Dependence of the half-life, r, on system size. v\\/v±_ = 1.764, as expected for 
DP in 2+1 dimensions. Filled squares: p = 0.3928; □: p = 0.3927; x: p = 0.3926; o: 
p = 0.3925. 

Fig. 2. Survival probability, mean number of open sites, and mean-square spread for 
natural initial conditions, p = 0.392615. 

Fig. 3. Local slopes, 5 t (a), rj t (b), and z t (c), for natural initial conditions, o: 
p = 0.39261; □: 0.392615; x: 0.39262; +: 0.39263. 

Fig. 4. Two-point correlation function h(r) in natural absorbing configurations. •: 
(1,0) direction; x: (1,1) direction; o: (2,1) direction. 

Fig. 5. Examples of initial absorbing configurations investigated in this work. Left: 
upper, density = 1/5; mid, = 1/4, anisotropic; lower, = 0.363, (RSA); right: 
upper, = 0.328, (natural); mid, = 0.25, isotropic; lower, = 0.5. 

Fig. 6. Local slopes, — 8 t (a), r] t (b), and z t (c), for = 0.2. +: p = 0.4106; □: 
p = 0.4104; o: p = 0.4102; x: p = 0.4100; o: p = 0.4096. 

Fig. 7. Spreading exponents 8, i], z, and (5' vs initial density 0. 

Fig. 8. Spread of open sites, = 0.2, p = p c . Light gray: t = 200; dark gray: t = 500; 
black: t = 1000. 
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Fig. 9. Mean density of open sites vs distance from seed, for = 0.2, p = p c . From 
left to right, t = 200, 500, and 1000. 

Fig. 10. A typical colony for = 0.5, p — p c , t — 10 4 . The line represents 100 lattice 
spacings. 
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